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The study of PDE's arose in the 18th century in the context of the devel- 
opment of models in the physics of continuous media, see [7J. It all began 
in 1747 when d'Alembert, in a memoir presented to the Berlin Academy, 
introduced and analyzed the one dimensional wave equation u# = u xx as a 
model of a vibrating string, see [T], and [2j. The appropiatness of this model 
still is a topic of discussion; for recent work in this direction, see [3], [6j, and 

d'Alembert observed that a general solution to the wave equation is given 
by u(x,t) = F(x + t) + G{x — t), where F, G are arbitrary functions, and 
F(x + t),G(x — t) represent waves that move along the string to the left 
and right, respectively, with constant speed. The Cauchy, or initial value, 
problem for the one dimensional wave equation consists of finding u(x,t) 
in the upper half-space given its initial displacement u(x, 0) = <p(x) 
and velocity ut(x,0) = ip(x). In terms of d'Alembert's traveling waves, the 
solution assumes the form 

Z Z Jx-t 

This expression serves as the prototype for the solution to the initial value 
problem for the wave equation in higher, odd, dimensions. 

The wave equation has a rich history fraught with controversy, both math- 
ematical, see [19] . physical, see [27], and philosophical, see [28J. At the heart 
of the controversy lies the notion of a function, as reflected by the initial 
shape of the string. Contrary to d'Alembert's view that only <p's and ip's 
that could be described analytically should be allowed as initial data, Euler 
believed that any curve, no matter how irregular, should be included, see 
[11] . At the same time D. Bernoulli, based on physical experiments, conjec- 
tured that the shape of a vibrating string, i.e., u(x,t), could be described 
by means of a trigonometric series, see [5], an assertion that would be vali- 
dated by Fourier's work on the heat equation, see [13]. Euler also believed 
that partial differential calculus should be extended to deal with irregular 
functions, a notion that was formalized in the twentieth century with the 
advent of generalized functions, or distributions, see [18] . [23] . These two 
ideas: the use of Fourier transforms and their interpretation as tempered 
distributions, are the essential ingredients in this paper. 
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Various methods have been utilized since to solve the Cauchy problem for 
the wave equation in higher dimensions - notably the finite part of a diver- 
gent integral, see [2], the analytic continuation of an integral of fractional 
order, see [21], non-self-adjoint equations, see [8], and the Poisson spherical 
means, see [20], [H], and [16] - see [9] for a nice account of these results. 

And, of course, the Fourier transform. By means of the Fourier transform 
in the space variables, a linear PDE with constant coefficients is reduced to 
an ODE in the time variable, which, once solved, gives the solution to the 
PDE in question by means of the Fourier inversion formula - the convolution 
properties of the Fourier transform and the fact that the Fourier transform 
of a spherically symmetric, or radial, function is again radial, are key ingre- 
dients in this context. Of the three major examples of classical second-order 
PDE's this approach provides a satisfactory solution to the Laplace and heat 
equations, but, for the wave equation the solution lacks the clarity of the 
other cases when the dimension n is greater than 3, see [10], [25], and |26| . 

In this paper I will derive the classical solution to the wave equation in 
dimensions n greater than or equal to 2 via the Fourier analysis techniques 
outlined above, thus putting it on an equal footing with the other equations. 

Upon embarking on this endeavor it becomes quickly apparent that in- 
tegrals of a function on M. n that depends only on x n need to be evaluated 
efficiently. Specifically, the following formulas will be used throughout the 
paper without further ado; for similar results, see |16| and [22j . 

Suppose n > 3. If the function f depends on x n only, 

f f(x n )dx = u n -i [ P [ f(s)(p 2 -s 2 ) ( - n ~ 3)/2 dsdp, 

JB(0,R) JO J-p 

and 

[ f(x n ) da(x) = u; n _! R (* f(s) (R 2 - s 2 )^' 2 ds , 
JdB(0,R) J-R 

where R > 0, da is the element of surface area of the sphere dB(0,R), and 
u n -i is the surface area of the unit ball ofM. n ~ . 

To verify these formulas one makes use of the n-dimensional polar coor- 
dinates given by 

xi = p sinOi) • • • sin(99 n _ 3 ) sin(v? n _ 2 ) sin(^ n _i) , 
x 2 = p sin(^i) • • • sin((/? n _ 3 ) sin(^ n _ 2 ) cos(v? n _i) , 
x 3 = p sin(^i) • • • sin(<£ n _ 3 ) cos((/? n _ 2 ) , 

x n _i = p sinOi) cos(v9 2 ) , 
x n = p cos(<^i) , 

where < (fk < ir, 1 < k < n — 2, < <p n -± < 2tt, and < p. As 
is readily seen by induction on n, the Jacobian of this transformation is 
J = p n ~ x sin n ~ 2 (93i) • • • sin 2 ((/? n _ 3 ) sin(^ n _ 2 )- 



Expressing the first volume integral above in polar coordinates gives 

r rR i"K 

\ f(x n )dx = c n / p n ~ l \ f (p cos((pi) ) sm n ~ 2 ((fx) d(fi! dp, 

JB(Q,R) J J 

where c n is equal to 

pn rir pin 

■ ■■ / sin n " 3 (99 2 ) • • • sin 2 O n _ 3 ) sin(</? n _ 2 ) d(f 2 ■ ■ ■ d<p n - 2 d(p n -i . 
Jo Jo Jo 

Now, this expression is readily recognized as the surface area uo n -i of the 
unit sphere in l"" 1 , or, equivalently, as the angular integrals correspond- 
ing to the evaluation of the volume v n —\ of the unit ball in R n_1 in polar 
coordinates, and so, c n = (n — 1) t> n _i = 0J n -i- 

As for the inner integral above, note that it equals 



f'TT 

/ /(pcos(pi)) sin n ~ 2 ((^i)(i99i 
Jo 



= [ W /(pcos(^) ) (1 - cos V) ) (n " 3)/2 sin(^!) fa , 
Jo 

which, with the change of variables pcos(<^i) = s, becomes 

r f(s) (i - {s/pf^-^p- 1 ds = r m ( P 2 - s 2 Y n - 3)/2 ds . 
j—p j—p 



So, combining these observations it follows that 

i-R 



I f(x n ) dx = w n _i [ p [ f(s)(p 2 -s 2 ) (n 3)/2 ds dp, 
Jb(o,r) Jo J-p 

which is the first relation above. The second is obtained at once from this 
by invoking the following simple consequence of Fubini's theorem: 

-R 



g(x) dx = / g(x) da(x) dr . 
Jo JdB(0,r) 



B(0,R) JO JdB(0,r) 

Indeed, differentiating with respect to R the identity 

f f f(x n )da(x)dr = uj n -! [ p f f(s)(p 2 -s 2 ) ( - n ~ 3)/2 dsdp, 

Jo JdB(0,r) Jo J-p 

it follows that 

/ f(x n ) da(x) = Wn _! R f R f(s) (R 2 - s 2)("- 3 )/ 2 ds , 
JdB(0,R) J-R 

and the formula has been established. ■ 

With this technical tool under the belt, the time has come to implement 
the Fourier analysis approach. First, some notations. 

A multi-index a = (a±, . . . ,a n ) is an n-tuple of non-negative integers. 
The length \a\ of a is \a\ = a\ + ■ ■ ■ + a n , and, for x = (x±, . . . ,x n ) G W 1 , 



x a = x^ 1 ■ ■ ■ x" n , and the differential monomial D a is defined as 

(d \ °i / d \ a n 
sr) ■■■(&-) 

Let 5(R ra ) denote the Schwartz class of rapidly decreasing functions, i.e., 
those functions <p on lR n such that, together with their derivatives, decay 
faster than any power of x as \x\ — >• oo. More precisely, ip £ S(R n ) iff 
suPiEgR™ \ x< * D l3( p(x)\ < oo for all multi-indices a,/3. 

Recall that the Fourier transform of a function if € 5(M n ) is given by 

m= I <p(x)e- ix <dx, 



and that it also belongs to In particular, note that D a ip(£ i ) 

iH£ a ^5(£). Furthermore, the Fourier inversion formula holds and 



Also, if if * ip(x) = j Rn ip(x — y) ip{y) dy denotes the convolution of tp and 
ip, it readily follows that ip * ip G <S(M n ), 

* vrco = m m , and, v * ^0*0 = 7^ / ft® $(® s ' x ■ 

As for the Fourier transform of a radial, or spherically symmetric, function 
(p on W 1 , it is again radial. Indeed, given a rotation matrix M in R n , since 
det (M) = 1 and Mx ■ M£ = x ■ £, it follows that 

<p(M£)= f ip(x)e- ix - M ^dx= [ l p(Mx)e- iMx - m dx = fi(0. 
Now, if 
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V = — — t: + ' ' ' + 



dx\ dx\ 

denotes the Laplacian in W 1 , the wave equation in the upper half-space 
is given by 

= V 2 u, ieK n ,t>0. 

ot 

The Cauchy problem for this equation consists of finding u(x, t) which 
satisfies the wave equation, subject to the initial conditions 

u(x, 0) = ip(x) , and ut(x, 0) = ip(x) , xeK", 

where, partly because of its simplicity, we restrict ourselves to functions (p, ip 
in the Schwartz class <S(R n ). 

Suppose u solves this problem. Taking the Fourier transform in the space 
variables, considering t as a parameter, since 

*Wl(£>*) = l<k<n, 



it readily follows that V 2 u(£, t) = — \C\ 2 u(^,t), and so, u satisfies 
u tt {^t) + \t\ 2 u(t,t)=V : ?GR n ,t>0, 

subject to 

S(f,0)=£(O, and £(£,0) (GR". 

For each fixed £ € M n this resulting ODE in t is the simple harmonic 
oscillator equation with constant angular frequency and one has 

*) = £(0 cos(t|e|) + ^(0 , e G M", t > . 

At least formally, the inversion formula gives for (x,i) € RIL + , 



u(x,t) = -±-[ coB(t|e|) e** < + 7^ / to^fP 



In fact, as long as (p and tp are sufficiently smooth so that differentiation in 
x and t can be performed under the integral sign, u(x, t) will be a solution 
to the Cauchy problem in the usual sense. However, a more transparent 
expression that does not involve the Fourier transform of the initial data is 
desirable, and, since the first integral above can be formally obtained from 
the second by differentiating with respect to t, I will concentrate on the 
latter. 

Now, the hope to interpret this integral as the inverse Fourier transform 
of a convolution involving ip is dashed because sin(|£|t)/|£| decays so slowly 
at infinity that it fails to be the Fourier transform of an integrable, or square 
integrable, function; this shortcoming is more apparent for cos(i?|£|), which 
does not even tend to as |£| — > oo. 

The next best thing is then to interpret sin(|£|i)/|£| as the Fourier trans- 
form of a generalized function, or tempered distribution. This will be a sim- 
ple consequence of two identities which, although classical in appearance, 
are not familiar to many. The difference in the character of these identities 
when the dimension is odd on the one hand and even on the other is not 
unusual in n-dimensional Fourier analysis; see [23] for further illustrations 
of this phenomenon. 

Specifically, the identities in M. n are given by the following expressions. 

Assume that n is an odd integer greater than or equal to 3. Then, 

ICI \RdRJ \UnRJdB(0,R) 

where R > 0, ui n is the surface measure of the unit ball in M. n , and c~ ] 
(n-2)(n-4)---l. 

On the other hand, if n is an even integer greater than or equal to 2, 

sin (^D =d (±±\ ^ 2 (J_ f 1 e -i*< dx 

1^ n \RdRJ \v n J B{m (fla-laf) 1 /* 



where R > 0, d^ 1 = n(n — 2)(n — 4) • • • 2, and v n is the volume of the unit 
ball in W 1 . 

Suppose first that n is an odd integer greater than or equal to 3. One 
begins by computing Xb(o,r)(Oi the Fourier transform of the characteristic 
function of B(0,R). Since Xb(o,R) i s radial, so is Xb(o,R)i an dj consequently, 
Xb(o,r)(0 = Xb(o,r) ((0, • • • , 0, -|C|)) • Thus, the integral can be expressed 
as 

/ / e~ ix< da(x) dr = Xb(o,R) (f ) 

JO JdB(0,r) 

= Xb(cvr)((0,...,0,-|£|)) 



= f e iXn ^dx 

JB(0,R) 



= Wn _! r P r e^\ ( P 2 - s 2 ) {n - 3)/2 ds dp . 

JO J-p 



-p 

Hence, differentiating with respect to R, 

[ e~ ix < da{x) = vj^R /"* e is ^ (R 2 - s 2 ) {n ~ 3)/2 ds , 

JdB(0,R) J-R 
or, in other words, 

[ R e is ^ {R 2 - s 2 ) {n - 3)/2 ds = I e~ ix < da(x) . 

J-R Wn-lR JdB(0,R) 

Note that, since n is an odd integer, m = (n — 3)/2 is also an integer. 
Furthermore, since (-R 2 — s 2 ) vanishes at R and —12, differentiating with 
respect to R gives 

2mR f R e^{R 2 -s 2 ) m - 1 ds = (—^— f e^<da(x)), 

J-R OR \0J n ^iR Job(0,R) ' 

or 

f R e «ICI f R 2 _ g2) m-l dg = 1 / 1 _1 /■ ^ N 

V ; 2m \R dR)\u n -\R JdB(o,R) J 

Clearly, repeating this process m times yields 

f R e^ds = ^-(^) m (^— I e-<da{x)). 
J-r 2 m m\\RdRJ \u n -iR Jqb(o,r) J 



Therefore, since 



it follows that 

sm(R\Z\) 



[ c n\t\ dc - 1 ( c iR\Z\ c -iR\&\ _ o s HR|g|) 

J-r m { j iei ' 

Ifl **Gias) (^R L(o,R) e ~ iX< da{x) ^ 



where m = (n — 3)/2, and 

Wn 1 1 



" n w n _i 2 2 m m! 
It only remains to compute c n . Since 

-m-2 



l — [ da{x) = R r 



'dB(0,B) 



setting |£| = in the relation just derived it follows that 

/ 1 d \ m „ 9 

and so, differentiating m = (n — 3)/2 times obtains 

i? = (n-2)(n-4)---l-c n i?, 

and the conclusion follows in this case. 

As for even dimensions, the proof follows from the odd case by essentially 
Hadamard's method of descent, see |14| . Since n is an even integer greater 
than or equal to 2, n+ 1 is an odd integer greater than or equal to 3. Observe 
that, for \ G M n , £' = (f,0) G E n+1 , and 

sin(^l^l) S in(fl|g'|) , _ p ^ n 

— iii — = — Hi — ' ^ eM ' jK>0 - 

Now, for x' G R n+1 , since x' = (x,x n +i) with x G M n , one has x'-£' = x-£. 
Moreover, in view that m' = (n + 1 — 3)/2 = (n — 2)/2, if B'(0,R) denotes 
the ball in M n+1 of radius i? centered at the origin, one has 

sin(i?|£|) _ f\ d \m' / 1 



K| c "+'(«m) fa/aw, 6 '"''J' 
Next note that dB'(0,R) can be decomposed into an upper and a lower 
hemisphere. The upper hemisphere, dB'(0, R) n {x^ +1 > 0}, is the graph of 

1 /2 

the function 7(x) = (i? 2 — |x| 2 ) for x G B(0,R), and, similarly, the lower 
hemisphere, dB'(0, R) n {x^ +1 < 0}, is the graph of —7. Furthermore, since 

V 7 (x)= 1 -— T7 ^x, x G M ra , 

(# 2 - |X| 2 ) 1/2 



it follows that 



R 



da(x') = vl + |V7(x)| 2 dx = -7^ fix 

(i? 2 -|x| 2 ) 1/2 

in each hemisphere. Hence, combining these observations, one has 

sin ^i) _ d (la.) (n - 2)/2 r 1 / 1 e -*< ^ 

|e| " n \R OR) \v n J mR) {R 2 _ | X |2)V2 

where 



To evaluate d n , we refer to the relation 



C n +1 



u n 2 2 m m'\ 
established above. Since nv n = u n , 

C n +1 V n 1 1 



d n = 2- 



uj n +i n 2 m m'\ 
Finally, in view that m! = (n — 2)/2, 

2 m 'm'l = 2("" 2 )/ 2 ( (n - 2)/2)! = ^-^^^^ (n - 2)(n - 4) • • • 2 , 

and, consequently, d n n(n — 2)(n — 4) • • • 2 = 1. ■ 

Finally, all the pieces are in place to carry out the program. Suppose first 
that n is an odd integer greater than or equal to 3. Two things must be 
proved: First, that the derivative of the surface area measure carried on the 
sphere dB(0, R) centered at the origin of radius R in R n given by 

T -c f 1 d \ { - n ~ S)/2 ( 1 f da(x)) 
\R dRJ \ u) n R Jqb(o,r) ' 

is a (compactly supported) tempered distribution, and, second, that 
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Both of these assertions can be verified painlessly. Recall first the notion 
of tempered distribution. A sequence {(fk} C <S(M n ) is said to tend to in 
5(R") provided that, for all multi-indices a, f3, the sequence {x a D^ifk} con- 
verges uniformly to for every x as k — >■ oo, i.e., limfc-^oo sup x \x a DPtpk{x) \ = 
0. A tempered distribution T is a linear functional T : S(R n ) — > R that is 
continuous in the sense that T(ipk) ~> whenever (fk —> in 5(R n ). 

As for the Fourier transform of a tempered distribution T, it is defined as 
the tempered distribution T that verifies T((p) = T(<p), for every ip E <S(R n ); 
that T is a tempered distribution is a consequence of the fact that, if {<Pk} 
tends to in 5(R n ), {<^k} also tends to in <S(R n ). Moreover, T*ip is defined 
for ip G 5(R n ), say, (T*<pJ^(£) = T(^) <p(£), and the inversion formula holds 
in this case as well, see [23], |25j . 

Now, to verify that Tr is a tempered distribution observe that, for ip £ 
5(R n ), 

T R (<p) =Cn (~) ( "" 3)/2 (R n - 2 — [ y{Rx') da(x' 



Since 



l<£<(n-3)/2. 



is a finite linear combination of terms of the form x' a D l3 ip(Rx') with |a|, \j3\ < 
£ and coefficients that are polynomials in R, a straightforward application of 



the Lebesgue dominated convergence theorem (LDCT) gives that the differ- 
entiation may be performed under the integral sign. Furthermore, if 
tends to in S(R n ), lim fe ^ oo sup^./ \x' a tp(Rx') | — 0, and, consequently, 
again by the LDCT, Tr(<^) — > as k — >■ oo. Hence, Tr is a (compactly 
supported) tempered distribution for each R > 0. 

Next, there are a couple of ways of computing Tr. First, for ip E <S(R n ), 




Now, since e l:E is integrable on R n x dB(0, R), by Fubini's theorem 
it is legitimate to change the order of integration, and so, 

T R {tp) = c n (I A) (n " 3)/2 f ( J_ f e -ix-€ d(7 ( x) ) v(0 ^ . 

\K OK/ J K n V W n it Jqb(0,R) j 

A direct argument along the lines of the one used to verify that Tr is a 
tempered distribution gives that 

JR n Is I 

and, consequently, Tr = sin(.R|£|)/|£| in the sense of distributions. 

On the other hand, one may note that, since Tr has compact support, by 
elementary properties of the Fourier transform of tempered distributions, 
Tr is the C°°(R n ) function given by Tr(£) = T R (e~ lx '^), and, consequently, 
again by the identity involving sin(i?|£|)/|£|, Tr(0 = sin(i?|£|)/|£|. 

Therefore, the integral of interest to us can be written as 

and so, by the inversion formula for tempered distributions, this expression 
is equal to T t *ip(x). So, it only remains to identify T t * ip(x). 

Since by definition T*7p is the C°°(R n ) function T*^(x) = t x T$), where 
t x denotes translation by x and ip(y) = ip(—y), it readily follows that 

/ 1 d \ (n-3)/2 /\f \ 

T t *il>(x) = c n (--) — - / ^{y)da(y)). 

\t dtJ \ U n tJdB{x,t) ' 

In other words, if 

M R <p(x) = l f i>{y) da(y) 
U n R n JdB(x,R) 

denotes the average of ip over the sphere centered at x of radius R, 

/ 1 8 \ (n-3)/2 



So, the solution to the Cauchy problem in this case can be expressed 
compactly in terms of spherical means as follows, see [12j . 

Let n be an odd integer greater than or equal to 3, and suppose that 99, ip 
are smooth functions. Then, with c" 1 = (n — 2)(n — 4) • • • 1, 

/ 1 d \ (n-3)/2 / „ 

is a C 2 (M™ +1 ) function that satisfies uu(x,t) = V 2 u(x,t), x G M n , t > 0, 
and u(x, 0) = nt(x, 0) = ip(x), x G M n . 

In /act, u;/ien u is defined as above in terms of spherical means, it suffices 
to assume that ip G C^ n+3 ^ 2 {R n ) and $ G C^ n+1 ^ 2 {R n ). 

When n = 3 the expression reduces to the familiar Kirchhoff 's formula. 

The result for even dimensions follows mutatis mutandis from the odd 
dimensional case. Indeed, if n is an even integer greater than or equal to 2 
and d~ l = n (n — 2) ■ ■ ■ 2, let 

r t = ,„(I*) ( "- 2,/2 (l/ 1 * 

\t at/ Vu n y B(ot) ^ 2 _ ixp) 1 / 2 

Then, T is a compactly supported tempered distribution, and Ti(£) = 
sin(t|£|)/|£|. Thus, if 

1 f 1 

v n R n Jb(x,r) [B 2 -\x- y\ 2 ) 1 
denotes the weighted average of ip over the ball centered at x of radius R, 

/I /J\(n-2)/2 

T t *^(x) = d n (-^) (t n M t if>(x)). 
Therefore, the following holds, see |12j . 

Let n &e an even integer greater than or equal to 2, and suppose that tp,ip 
are smooth functions. Then, with d~ l = n (n — 2) • • • 2, 

d /l S\M/a, 



1 (9 \ (n-2)/2 

+ ^(7^7) (t n W(*: 



.t dt, 

is a C 2 (M™ +1 ) function that satisfies u tt {x,t) = V 2 u(x,t), x G M n , i > 0, 
and n(x, 0) = <p(x), ut(x, 0) = ip(x), x G R". 

In fact, when u is defined as above in terms of means, it suffices to assume 
that v G C( n+4 )/ 2 (M n ) and ip G C( n+2 )/ 2 (M n ). 

In particular, when n = 2, the result is known as Poisson's formula. 

I am grateful to the referee and the editor for their valuable comments 
that shaped the final presentation of this paper. 
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